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ABSTRACT 

Modeling the multiwavelength emission of successive regions in the jet of 
the quasar PKS 1136-135 we find indication that the jet suffers deceleration 
near its end on a (deprojected) scale of ~400 kpc. Adopting a continuous flow 
approximation we discuss the possibility that the inferred deceleration from a 
Lorentz factor T = 6.5 to Y — 2.5 is induced by entrainment of external gas. 

Some consequences of this scenario are discussed. 
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1. Introduction 

Despite decades of intense efforts, the present knowledge of the physical processes acting 
behind the phenomenology shown by relativistic jets is still rather poor. Basic questions 
about matter content, transported power, the role of the magnetic field, the dissipation 
mechanisms are awaiting answers (e.g. Blandford 2001). Among these problems, one of the 
most fundamental concerns the speed of the flow. While it is widely accepted that the rapid 
variability and the extreme compactness displayed by blazars (e.g. Maraschi & Tavecchio 
2005; Sikora & Madejski 2001) and the observation of VLBI superluminal components (e.g. 
Kellermann et al. 2004 ) imply relativistic bulk flows (with T ~ 10 — 20) for the innermost 
portions of the jet (0.1 pc < d < 100 pc), it is not clear whether the large scale jets still have 
relativistic speeds and how/where deceleration takes place. The present evidence suggests 
that FRI jets become trans-rclativistic quite early, within few kiloparsecs (e.g. Bridle & 
Perley 1984), while the situation for FRII jets appears more ambiguous. Mildly relativistic 
speeds (/? < 0.95) are suggested by studies of high luminosity, lobe dominated sources 
based on the jet to lobe prominence (Wardle & Aaron 1997), while the interpretation of 
multiwavelength observations of extended jets in QSOs points toward highly relativistic 
speeds (T ~ 10) even at very large scales (~ 100 kpc; e.g., Tavecchio et al. 2000, Celotti et 
al. 2001, Sambruna et al. 2002,2004, Siemiginowska et al. 2002, Marshall et al. 2005). 

Deceleration of FRI jets is commonly believed to be the result of entrainment of external 
gas by the flow (e.g. Komissarov 1994, Bickncll 1994, B94 hereafter). In this framework it is 
conceivable that FRII jets, characterized by larger powers, are not substantially perturbed 
by entrainment: the fundamental division between FRI and FRII jets could basically depend 
on the difference in the mass/energy flux of the two types of jets (Bicknell 1995, Bowman 
et al. 1996). An important unknown parameter is the efficiency with which the jet can 
entrain the external gas. Efficient entrainment can be assured by the onset of fully-developed 
turbulence at transonic speeds (e.g. B94). In this scheme, supersonic flows are thought to 
be less efficient at collecting material and entrainment could affect only the external regions 
of the jet, resulting in the formation of a slow layer surrounding a fast, unperturbed “spine” 
(e.g. de Young 1986). Nevertheless, even low efficiency entrainment could become important 
if it is maintained over very long scales. 

The possibility of “observing” the gradual slowing-down of a jet could in principle pro¬ 
vide precious information on the physical processes at work. This approach was successfully 
developed in great detail for a few FRI jets where the morphology could be well studied 
thanks to the large angular scale (Laing & Bridle 2002, 2004). For FR II jets the stronger 
asymmetry and smaller angular scale has generally prevented comparable studies. 

Of particular interest are therefore the cases of FRII jets hosted by intermediate-z 
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QSOs in which multiwavelength observations have recently shown that the radio to X-ray 
flux ratio increases dramatically along the jet. This feature, interpreted in the framework 
of the synchrotron-IC/CMB emission model, suggests that the jet undergoes deceleration. 
The possible role played by deceleration was first discussed for the jet of 3C273 (Sambruna 
et al. 2001). Other cases of jets showing an increasing radio to X-ray flux ratio were 
presented in the survey discussed in Sambruna et al. (2002) and Sambruna et al. (2004). 
Georganopoulos & Kazanas (2004; hereafter GK2004) proposed an analytic framework to 
describe the effects of deceleration on the physical parameters of the jet on the basis of an 
adiabatic assumption, remarking that, if the jet decelerates, the magnetic field as well as the 
particle energy density will increase due to compression causing a large increase of the radio 
synchrotron emission, while the X-ray flux will decrease because of a reduced beaming of 
the CMB radiation field. In Sambruna et al. (2005; Paper 1 in the following) we report the 
results of deep Chandra and multiwavclenth observations of PKS 1136-135, which allowed 
us to image in detail the emission along the jet. Modelling the observed multifrequency 
emission we derived the profiles of the basic physical quantities along the jet, which support 
the deceleration scenario discussed above. In this paper on the basis of the observational 
evidence collected for the jet of 1136-135, already presented in Paper I, we discuss in more 
depth deceleration in terms of physical models and in particular entrainment of external gas. 

The plan of the paper is the following: in Sect.2 we derive the profile of the relevant 
physical quantities, suggesting deceleration of the flow. In Sect.3 we discuss the physical 
origin of the deceleration. In Sect.4 we discuss in detail the possibility that deceleration 
is induced by the loading of the jet through entrainment. Finally in Sect.5 we conclude. 
Throughout this work we use the current “concordance” cosmological parameters: Ho = 
71 Km s^ 1 Mpc -1 , Ha = 0.73, H m = 0.27 (Bennett et al. 2003). 


2. Determining the deceleration profile 

One of the most striking features of the the jet of 1136-135 after knot B is the dramatic 
increase of its radio luminosity and the associated decay of its X-ray emission (Fig.l). The 
brightness profiles directly translate into the run of the physical quantities found in Paper 
I by using the IC/CMB radiative model (Tavecchio et al. 2000; Celotti et al. 2001). In¬ 
deed the increasing radio brightness requires an increase of one order of magnitude both of 
the magnetic field and the density of the radiating electrons (bound together through the 
equipartition condition). The decrease shown by the X-ray flux, in turn, implies a decrease 
of the Doppler factor, which reduces the amplification suffered by the target photons of the 
CMB. Therefore, the observed morphology and the resulting profiles suggest that the special 



4 


behavior displayed by 1136-135 could be related to effects induced by the deceleration of the 
flow, starting in the region corresponding to knot B-C. Remarkably, the physical parameters 
inferred for the hot spot (F in the prsent work, HS in Paper I) appear to continue the same 
trend followed by the knots, supporting the emerging view that the hot spots of jets in pow¬ 
erful QSOs can be thought of as slow knots, rather than portions of the jet almost at rest 
(Georganopoulos & Kazanas 2003; Tavecchio et al. 2005). We remark that the morphology 
and the results of the modelling suggest that the deceleration is important only after the 
region B-C. Before that point the character of the flow seems to be different, with emission 
clearly concentrated in knots characterized by a similar speed. In this paper we only consider 
the region after knot B, where deceleration seems to take place. 

We review the main features of the model used to fit the data. The emitting region is 
assumed to be in motion with a bulk Lorentz factor T at an angle 6 with respect to the line 
of sight. Since fluxes are extracted within elliptical regions (Paper I), we assume a volume of 
the emitting region equal to that of the corresponding (deprojected, see below) ellipsoid. The 
emitting region is homogeneously filled by high-energy electrons, with a power-law energy 
distribution iV(y) = K r )~ n (K is the electron normalization) extending from 7 m i n to 7 max . 
Electrons emit radiation through the synchrotron and IC/CMB mechanisms. The low-energy 
limit of the electron distribution 7 min is constrained by the optical and the X-ray fluxes. The 
calculations reported in Paper I have been performed assuming that the observed emission 
concentrations (knots) mark the presence of discrete moving components of the jet. In the 
region of the jet under consideration here (after knot B) the flow seems to become more 
similar to a continuous flow, in which it is difficult to isolate knots (Paper I). It is likely 
that the flow is complex, showing the characteristics of a continuous flow together with 
the presence of regions where the dissipation is more efficient, visible as brighter regions in 
the jet. Based on this evidence, and to be consistent with the scenario that we intend to 
explore in this paper, we recalculated the physical quantities assuming that, after knot B, 
the emitting plasma behaves as a continuous flow, moving in a jet whose pattern is at rest 
with respect to the observer. Practically this implies that the exponent 2 + a should be used 
instead of 3 + a for the Doppler factor to describe the amplification of the radio flux (Lind 
& Blandford 1985) and that the observed length is deprojected using the observing angle. 
The choice of this particular geometry slightly affects the values of the derived quantities, 
but does not change the evidence of a trend suggestive of deceleration. 

With the assumptions above it is possible to derive the required values of the physical 
quantities by using the observed radio (5 GHZ) and X-ray (1 keV) fluxes. To break the 
degeneracy between the model parameters (we have three physical quantities the magnetic 
held, B, the electron normalization K and the Doppler factor 5 and two quantities from ob¬ 
servations, the radio and the X-ray flux) it is necessary to introduce an additional condition. 



A widely adopted assumption (although, admittedly, not fully motivated) is that particles 
and magnetic field are in equipartition (e.g. Tavecchio et al. 2004; but see Kataoka and 
Stawarz 2005). Several evidences indicate the presence of a proton component in jets (e.g. 
Sikora & Madejski 2000). However, we assume that protons are not directly coupled to the 
non-thermal components (relativistic electrons and magnetic field), although they can con¬ 
tribute to the total pressure. We then assume equipartition between relativistic electrons and 
magnetic field energy densities, Ub — U e , assuming no contribution from protons. Note that, 
in any case, the inclusion of protons into the equipartition condition only slightly affects the 
derived quantities (unless protons strongly dominate the total pressure). It is then possible 
to find analytical relations providing a unique set of physical parameters characterizing the 
emission region (e.g. Tavecchio et al. 2000). The appropriate analytic relations are given in 
Appendix A. 

We devoted special attention to the estimate of the uncertainties affecting the final 
parameters (B, K and 5) deriving from the uncertainties associated to the input parameters 
of the model. To this aim we calculated the ranges covered by B, K and 5 for different 
values of the input parameters, varying j m i n in the range 5-25, the volume of the emitting 
region in the range V — V/2, the radio and the X-ray fluxes and the radio spectral slope 
a r within the corresponding errors (given in Paper I). We find that the total uncertainty 
is largely dominated by the error on the radio spectral index. We then associate the total 
range spanned by B } K and 5 for the different input parameters with the total uncertainty. 

A more complex discussion involves the viewing angle, 6. It enters as an input parameter 
in the calculations, through the deprojection of the observed length. At the same time, the 
angle is constrained by the results of the calculations. Indeed, for a derived Doppler factor 

5 (resulting from the modelling), 9 cannot be larger than 1/5 (see e.g. the discussion in 
Tavecchio et al. 2004). To fix the angle we then follow an iterative procedure: we fix an 
angle, we derive the parameters, we check if the maximum Doppler factor required by the 
knots (including the upper bound deriving from the uncertainties discussed above) satisfies 

6 < 1 /5 max . If the condition is not satisfied we reduce the input viewing angle and we repeat 
the procedure, until we reach the maximum angle allowed by the data, 0 max = 3.8 deg. All 
the angles 6 < 9 max are in principle allowed. However, since # max is the most probable angle 
and the impact of a different choice is quite small and 9 max is the most probable angle, we 
use # max in the modelling. 

The profiles of B , K and T, derived from the emission model with equipartition, are 
shown in Fig. 2. The value of the Lorentz factor T can be derived from 5 once the observing 
angle is fixed. We recall that each knot is modelled separately from the others: therefore 
the trends we derive are not by-products of the procedure used, but characteristics of the 
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flow, and are basically related to the observed trends of radio and X-ray fluxes along the 
jet. The evolution of the model parameters with distance are in qualitative agreement with 
the adiabatic model discussed in GK2004. Unfortunately, the large uncertainties prevent a 
deeper discussion of this point. 

In conclusion our results support the view that the jet of 1136-135 suffers strong decel¬ 
eration after knot B. The Lorentz factor, T changes, from (with this specific choice of the 
angle) 6.5 to 2.4 in a projected angular distance of ~ 4 arcsec, corresponding to a deprojected 
linear distance of ~ 400 kpc. 


3. Origin of deceleration 

In the following we consider possible deceleration mechanisms, a question not explicitly 
addressed in GK2004. First we estimate the possible role of the Compton-Drag effect and 
then proceed to discuss entrainment. 


3.1. Compton Drag 

A possible way of decelerating a relativistic plasma is through the recoil induced by 
Compton scattering between the electrons present in the jet and some external radiation 
held (a process known as Compton Drag effect, e.g. Sikora et al. 1996). This mechanism 
has been recently proposed as a viable way to decelerate BL Lac jets within the first parsec 
from the core (Ghisellini, Tavecchio & Chiaberge 2005) and as a source of IR radiation from 
interknot regions in large scale jets (Georganopoulos et al. 2005). 

In the particular case analysed here the external radiation held is dominated by the 
CMB (with energy density Uqmb), while contribution of radio emission from the slow part 
of the jet (beamed in fast jet frame) is at most comparable to that of the CMB (e.g. Celotti 

et al. 2001). A key role in determining the efficiency of such a process is played by the 

composition of the jet: clearly, only when electrons carry a large fraction of the total jet power 
the process can affect the jet dynamics. Basically, the equation describing the evolution of 
the bulk Lorentz factor with time can be written as (e.g. Ghisellini et al. 2005): 

• <jtc < q 2 > T 2 £/cmb 

1 -7 w l 1 ) 

2 (u^+ < 7 >) 

where u = n p /n e . < q 2 > and < 7 > depend on the detail of the electron distribution. 
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Using typical values derived with the IC/CMB model, n = 2.6, 7 m i n = 10, 7 max = 10 5 (e.g. 
Sambruna et al. 2004), they can be estimated < y 2 >~ 10 4 and < 7 >~ 25. We assume the 
most optimistic case in which electrons do not substantially cool (or are continuously heated) 
and therefore < y 2 > is constant. The effects of the deceleration can be estimated through 
the quantity r ( \ ec = cT/T, yielding the typical distance at which the deceleration takes place. 
Inserting the numerical values we find rd ec ~ 100 Mpc for proton-dominated jets [u — 1), 
t ,dec ~ 300 kpc for pair-dominated jets (u> = 0) with < 7 >~ 10. We conclude that Compton 
Drag can be efficient only in the case of a pair-dominated jet. General arguments applied to 
blazars (although not completely conclusive) seem to point toward a proton component for 
jets in quasars dominating the dynamics (e.g. Sikora & Madejski 2000, Ghisellini & Celotti 
2002a, Maraschi & Tavecchio 2003; see also Uchiyama et al. 2005). We can conclude that an 
important role of the radiative recoil in determining the inferred deceleration seems unlikely, 
although not completely ruled-out. 


3.2. Entrainment 

In the following we discuss the possibility that the inferred decline of T along the jet 
is caused by the continuous matter loading of the jet through entrainment of external gas 
entering from the lateral walls of the jet. This mechanism is believed to be responsible for 
the early deceleration of FR I jets (e.g. B94). 

Basically, deceleration through entrainment can be understood to happen through a 
continuous series of inelastic collisions between the moving plasma and the external gas at 
rest (e.g., Icke 1991). As a result of the collision, part of the kinetic energy is dissipated 
and converted into internal energy of the jet, thus increasing the internal pressure. We first 
derive an approximate formula for the evolution of the Lorentz factor in the simplified case 
of an entraining confined mass (cylinder). Although drastically simplified, this approach 
allows to identify the key elements of the entrainment process. 

During the motion, the external layer of the cylinder interacts with the external medium, 
resulting in a continuous increase of the moving mass. Hence a continuous deceleration clue 
to momentum conservation and an increase of internal energy clue the dissipation of kinetic 
energy follow. The dynamics of the entraining mass can be described using the conservation 
laws for energy and momentum for inelastic collisions. The initial mass is M a with initial 
Lorentz factor is T 0 . Assuming that the entrained gas is initially at rest, conservation of 
energy and momentum yields: 
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where we have indicated with m the mass of the material entrained by the blob and E is 
the internal energy accumulated in the collision. Eqs.(2)-(3) give the Lorentz factor T as a 
function of the rest-mass m of the collected material: 


r 0 M 0 + m = YM = Y [ m + 


r 


a + r 0 

(a 2 + 2 aT 0 + l) 1 / 2 


( 4 ) 


where a = m/M 0 . The Lorentz factor reached depends only on the ratio between the total 
entrained and initial mass of the blob. For m/M 0 « l/r o the Lorentz factor remains almost 
constant, while for m/M 0 » 1/Y 0 Y tends to 1. Thus the critical scale for deceleration is 
reached when the collected mass is of order m crit ~ M 0 /Y 0 . To halve the initial bulk Lorentz 
factor T 0 the collected mass must be m = 1.5M 0 /T o . In the limit of the inelastic assumption 
the collected mass needed to slow down a highly relativistic mass is smaller for higher T 0 . 

In the following we apply the rigorous hydro dynamical treatment developed by B94 
to describe the deceleration of the jet of 1136-135 in order to discuss the plausibility of 
the entrainment mechanism for this particular case. Laing & Bridle (2002) used the same 
approach to model the dynamics of the jet of the FRI galaxy 3C31. Our analysis will be not 
as detailed, due to the more limited information allowed by the larger distance and smaller 
angular scale. The treatment assumes that the system is adiabatic, the dissipated kinetic 
energy being stored in the plasma internal pressure. This condition is satisfied in the case 
under study, for which it is known that the radiative efficiency is quite low (e.g. Tavecchio 
et al. 2000; Sambruna et al. 2002). 

We use momentum and energy conservation in the form given by B94 to follow the 
evolution of the decelerating flow. An important parameter is y = pc 2 /(e + p), representing 
the ratio between the rest mass density of the particles within the jet and the sum of internal 
energy density and pressure (enthalpy). Basically this parameter measures the ratio between 
the rest-mass energy and the (random) kinetic energy of the particles of the jet, as measured 
in the comoving frame of the flowing plasma. Values of y > 1 indicate cases for which 
the average kinetic energy of the particles is less than the rest mass-energy, i.e. a non 
relativistic plasma. On the contrary values of y < 1 characterize a relativistic plasma. 
Note that, based on the modelling of the observed emission, we can directly probe only the 
non-thermal component of the plasma (relativistic electrons and magnetic held). However, 
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the plasma can also contain non relativistic or mildly relativistic electrons and protons (the 
“thermal” component), contributing to (and possibly dominating) the mass density and the 
total pressure. Defining the initial internal pressure of the jet p 0 > the initial Lorentz factor r 0 , 
the initial cross sectional area of the jet A 0 , the initial y 0 parameter and the corresponding 
downstream parameters, the coupled expressions for momentum and energy conservation 
read: 


and: 


(i + X )r 2 p 2 + \ 


(1 + Xo)To/3o + 



[(r - i) x + r] r/3 = [(r 0 - i)y 0 + r 0 ] r 0 p 0 



(5) 

( 6 ) 


B94 assumed that the dominant role in determining the pressure is played by relativistic 
particles for which the relativistic equation of state p = e/3 can be applied; we adopt the 
same assumption. As can be verified a posteriori , the assumption is fully consistent with 
the results only at the end of the deceleration process. However, we verified that even in 
the unrealistic case of a non-relativistic equation of state there are just small corrections to 
the results (see also Laing & Bridle 2002) and the general conclusions of our study do not 
change. Analogously, for simplicity, in Eq.(5) the terms corresponding to the effects of the 
external medium (expected to be small) have been dropped. 

To solve the equation above we assume that the cross sectional area of the jet, A, remains 
approximative^ constant from C to F. Although this assumption is not fully consistent with 
our choice to neglect the effects of the external medium, the “cilindrical” approximation 
seems tolerable since the expansion factor in the considered region is at most a factor of 2. 
Moreover the jet could be self-confined by, for instance, a large scale magnetic field. 

Once the initial state is defined, the coupled non-linear equations can be solved to 
provide p and x a t each value of T, considered as a parameter. Other important derived 
quantities, such as the comoving density p and the entrainment rate, can also be calculated. 


4. Results 

The basic features of the deceleration process are shown in Fig. 3 through the evolution 
of the parameter y, the density and the pressure as a function of the Lorentz factor reached 
by the flow. For this calculation we assume initial conditions corresponding to the parameters 
found above for knot C of 1136-135, after which the deceleration is more evident (Fig. 2). 



As discussed above, we do not have any information about the thermal components of the 
plasma. The mass density has been calculated assuming a composition of 1 cold proton per 
emitting electron. For the pressure at knot C, we assume two cases. In the first, we assume 
that the pressure is dominated by the non-thermal components, p Q = p e +PB (solid lines). 
In the second case (dashed lines) we assume a pressure ten times that the non-thermal one, 
thought to be provided by a thermal component. 

To satisfy energy and momentum conservation, each small entrained mass induces a 
small deceleration of the flow and a small increase of the internal energy (assumed to be 
shared among all the protons inside the jet). This produces variations of the y parameter, 
as shown in Fig. 3 (top). Assuming a jet initially dominated by cold gas y >> 1, the 
dissipation induced by the continuous loading implies a decrease of y, because the pressure 
increases faster than the density of the jet. When the jet is almost completely stopped y 
is of order unity and starts to increase. Note that the regime considered by Laing & Bridle 
(2002) for 3C31 corresponds to low T, corresponding to the region where y increases during 
the deceleration. In our application to a highly relativistic FRII jet, we are considering the 
regime where y is large at the beginning and decreases along the jet. 

The profile of the comoving density (normalized to the initial value) is shown in Fig.3 
(middle). Also shown in the same figure is the profile one would get in the case of deceleration 
without entrainment (with the same geometry used for the other calculations): in this case 
the increase of the density only reflects the effect of deceleration and consequent compression 
of the flow (following the law p/po = r 0 /r expressing the conservation of the number of 
particles). The difference between this curve and the other curves represents the contribution 
to the total density of the entrained material accumulated along the path. Comparing the 
densities with and without entrainment it can be seen from Fig. 3 that the entrained mass 
(proportional to the density difference) needed to decelerate the jet from r 0 ~ 6.5 to V ~ 3.2 
(halving the initial Lorentz factor) represents a small fraction (about 10%) of the original 
rest-mass. 

In the bottom panel of Fig.3 we report the pressure of the plasma for different Lorentz 
factors reached by the flow. We also plot (points) the values of the pressure supported by 
the non-thermal components only (magnetic field and non-thermal electrons), as estimated 
from the modelling of the emission. The errorbars along the a:—axis indicate the uncertainty 
associated to the Lorentz factor of each region, as reported in Fig.2. Clearly the non-thermal 
component alone can account for a small fraction (around 10%) of the pressure required by 
the conservation laws, even in the case in which the initial pressure is completely non-thermal. 
Therefore much of the pressure is provided by the energy stored within the invisible thermal 
component. 
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The ability of the entrainment process to decelerate the jet is related to the quantity 
of matter effectively entrained by the jet along its path, expressed by the entrainment rate. 
The essential parameter characterizing this quantity is the entrainment velocity. Attributing 
the “observed” deceleration to entrainment we can derive the required entrainment rate and 
entrainment velocity in the specific case of 1136-135. 

The entrainment rate can be defined starting from the law describing the evolution of 
the jet mass flux (see Eq. 9 of B94): 

r 2 p 2 /3 2 cA 2 = T 1 p 1 (3 1 cA 1 + p ext VentdA (7) 

where subscripts 1 and 2 indicate quantities of the jet measured at distances x\ and x 2 , 
dA ~ 2nrjdx is the jet lateral surface and u e nt is cm effective entrainment speed. The 
(average) entrainment rate is thus defined as: 

r 2 p 2/ d 2 cA 2 — TipiPicAi If" 

- At - = At / P** V «* dA 

^v£l,2 cA-£i,2 


The average entrainment rate calculated from C to F (comparable values can obtained 
for the intermediate regions C-D, D-E, E-F) is ~ 2 x 10 24 g s~ 4 kpc -1 (the distance has 
been deprojected assuming 9 = 3.8 deg; the inferred entrainment rate would scale as oc 9). 
For comparison, this value is about an order of magnitude larger than the level required for 
the deceleration of the jet of 3C31 (Laing & Bridle 2002). In the latter case the collecting 
surface is smaller and the density much higher. Further assuming that the external density 
and the entrainment velocity do not change dramatically from C to F and specifying the 
external density, it is possible to derive the entrainment speed v ent . Assuming that in the 
region from C to F the external gas has a density n ex t = 10” 5 cm -3 (see below) and that the 
jet in the region considered here has a radius of r = 10 23 cm we find u ent = 4 x 10' cm/s. 
The corresponding value for the jet of 3C31 lies around u ent = 10 5 — 10 6 cm/s. Therefore, 
the entrainment process in 1136-135 appears to be more efficient than that occurring in the 
low-power jet of 3C31. 


4.1. The onset of deceleration 

A question naturally arising is whether the parameters required by the entrainment 
model for the ’’observed” deceleration in the region C-F are consistent with the propagation 
of the jet from the nucleus to the large scale. For this scope, we consider the simplest 
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hypothesis that entrainment is active with the same efficiency (measured by the parameter 
n en t) all along the path of the jet and assume a conical, entraining jet 1 . 

It is possible to check that the critical “cumulative” entrained mass flux for substantial 
deceleration in the case of a continuous jet satisfies a relation analogous to that derived 
above for the case of a blob, m CI it = M 0 /T 0 . Deceleration becomes important when the 
contribution of the cumulative mass entrained by a jet slice along its path is of order 1/T of 
the jet slice mass. This condition can be put on an explicit form: 


Pext^'enl-bl 


" AT,1 

tt 


= pifiicAi 


(9) 


In order to estimate the entrained mass flux (Eq.9) we need to specify the distribution 
of the external gas density, p ext — n ext m p while v ent will be assumed to be constant. From 
the Chandra image we infer the presence of faint diffuse emission surrounding the QSO. The 
distribution of the gas can be modelled as a King profile with core radius x Q = 40 kpc and 
central density n a — 5 x 10 -4 cm -3 (Sambruna et al., in prep). This result is consistent 
with other X-ray observations which show that, in general, nearby (z < 1) QSOs lie in 
environments typical of small groups (Hardcastle & Worrall 1999; Crawford & Fabian 2003). 
At distances of interest here the density can be well modelled as a power-law: 


n(x) = n Q 



( 10 ) 


with x 0 = 40 kpc, n 0 — 3 x 10 -4 cm -3 . From the profile above, the density in the region from 
knot C to F is expected to be of the order of n ext ~ 10~ 5 cm -3 , which has been used above 
to derive the value of u ent ~ 10 8 cm/s. Assuming that u ent does not substantially change 
along the jet, it is thus possible to estimate the total matter collected by the jet from the 
kpc scale to the hundreds of kpc scale, where deceleration becomes effective. Assuming a 
conical jet, the area dA can be expressed as dA = 2nrdx, with r = r 0 x/x 0 , where r Q is the 
radius at the initial point x Q (typically it is assumed r 0 /x 0 ~ 0.1). 

Even if the entrainment efficiency is constant and the density of the external medium is 
decreasing, the entrained mass flux steadily increases, due to the continuous increase of the 
jet’s surface. Calculating the integral in Eq. (9) we ford, as expected, that the typical scale 


Hn the previous section the jet radius has been approximated as constant in the region from C-F. However, 
in the following we want to study the cumulative effect of entrainment from very small to large scales, and 
thus the increase of the jet radius with the distance must be taken into account. 



13 


at which deceleration is effective is x\ ~ 100 kpc. Thus the propagation of the jet up to large 
distances can be understood. The terminal part of the jet of 1136-135 can be identified as 
the region where the entrained mass becomes critical causing substantial deceleration. We 
note however that in the schematic model we have drawn with constant entrainment speed, 
the deceleration profile is very regular. Clearly there must be additional conditions (e.g. 
the onset of strong instabilities) which enhance the process at particular sites driving the 
localized occurrence of deceleration. 


5. Discussion and Conclusions 

The proposed scenario can explain in a plausible way the deceleration inferred for the jet 
of 1136-135 (and possibly in the other sources showing the same radio-to-X-rays increasing 
trend, GK2004). A more detailed understanding of the processes at work and the comparison 
with the observed properties of jets necessarily involves several still poorly-known physical 
issues. 

It is worth remarking that the evidence for deceleration in the jet of 1136-135 is based 
on the observed run of radio and X-ray fluxes along the jet and on the use of a specific 
emission model to interpret the multiwavelength observations (for a discussion of alterna¬ 
tive emission models see Paper I). In this respect, the conclusion about deceleration is not 
completely model-independent. Other specific assumptions made in the modelling, such as 
the equipartion condition, affect (although in a minor way) the conclusion. Moreover, apart 
from the application of a specific emission model, one has to bear in mind that the result is 
based on the hypothesis that the emission comes from a not-structured flow. At least part 
of the increase of the radio flux close to the terminal part of the jet could be due to the 
contribution from another component (e.g. the plasma backflowing after the compression at 
the hot-spot). 

The approach that we used to study the deceleration is based on a pure hydrodynamical 
treatment. The contribution of all the components of the flow (protons, electrons, magnetic 
field) is included in just one parameter, the total pressure. Likely, much of the energy dis¬ 
sipated during the deceleration will be stored in the heated proton component. To follow 
the evolution of the other components (thermal and non-thermal electrons, magnetic field) 
we should specify the actual mechanisms coupling each component to the others. The de¬ 
termination of the coupling between protons, magnetic field and electrons, determining the 
fraction of dissipated energy that goes into the magnetic and non-thermal electron compo¬ 
nents, is a longstanding problem in astrophysics, clearly beyond the scope of this work. The 
derived profiles are in agreement with the possibility that the increase of the magnetic held 
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and of the density of the relativistic electrons just follows the adiabatic heating induced 
by the deceleration. However, given the large uncertainties, this evidence is not conclusive. 
The contribution of non-thermal electrons and magnetic held to the total pressure is small, 
around 10%. Therefore, a large portion of the dissipated kinetic energy resides in the hot 
protons and is not radiated away. 

Another interesting point is related to the nature of the mechanism able to induce the 
assumed entrainment. It is widely assumed that turbulence can effectively work only in 
slow (transonic or subsonic) jets. Apart from fully developed turbulence, it is also possible 
that entrainment is mediated by viscosity effects induced by small-scale turbulence. It is 
worth noting that in the calculation we have presented it is supposed that the entire jet 
is decelerated under the effect of the entrainment. More realistically, these processes will 
induce a radial structure in the jet, with the development of a sheared how (e.g. De Young 
1996). 

In the simple scheme presented here, the jet starts to be decelerated when the amount of 
collected gas is of the order of 1/T of the mass transported by the jet. In this framework, jets 
characterized by different mass huxes will experience different behaviours. Large mass fluxes 
will assure that, under the same conditions of external gas density and entrainment rate, the 
jet will reach its hotspot almost unperturbed. On the other hand, jets with a small mass hux 
will be decelerated soon. It is tempting to further speculate along these lines, connecting the 
deceleration of the jet to other properties of the central AGN. To this aim we indicate the 
jet mass flux as M out , and we note that M out is the mass output rate of the central engine, 
which is probably some fraction of the accretion mass rate M out oc M acc (e.g. Ghisellini & 
Celotti 2002b). From the condition expressed by Eq.(9) and the definitions above we can 
find that the deceleration length x\ is x\ oc M acc /r 0 , (the exact dependence is related to 
the jet and the external density profile) where we have supposed that all the jet, before the 
deceleration, have a similar bulk Lorentz factor. If the initial radius of the jet scales as the 
gravitational radius of the central black hole r Q oc Mbh, we finally find x\ oc m acc , where 
m acc is the accretion rate in Eddington units. Thus it is possible to identify sources for which 
the jet is effectively decelerated within short scales with sources with low m acc and viceversa, 
a situation reminiscent of the FRI/FRII division. 


We thank Gabriele Ghisellini for useful discussions and the anonymous referee for con¬ 
structive comments and criticisms. F.T and L.M acknowledge support from grant COFIN- 
2004-023189-005. 
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A. Analytical estimate of the physical parameters 

• Synchrotron EC/CMB emission. Assuming a simple power-law distribution for the 
emitting electrons: 

A( 7 ) = Kr n , (Al) 

the synchrotron luminosity can be written as (e.g. Ghisellini et al. 1985): 


L s {y s ) = c{o)KB l+a V8 p v~ a (A2) 

where V is the volume of the emitting region, a — (n — 1 / 2 ) is the spectral index, c(a) is a 
constant and the luminosity is calculated at the frequency (observer frame) v s . The exponent 
p of the Doppler factor depends on the geometry of the flow: p = 3 + a for a moving “blob”, 
p — 2 + a for a flow whose pattern is fixed in the observer frame (Lind & Blandford 1985). 

To express the luminosity emitted through the IC/CMB ee start from the the standard 
relation for the inverse Compton luminosity: 


A Hmax 

Ldy) = t^vtcVU , ad £ Pl / 7 2 N(p)5(u - v c )dp (A3) 

^ 'Tmin 

where pi = 4 + 2a for a “blob” (Dernier 1995) and pi = 3 + 2a for a continuous flow. The 
above expression assumes that the spectrum of the soft radiation is extremely peaked around 
the frequency u Q and then 7 and v c can be related by: 



Using this expression we can rewrite eq. (A3) as: 

L c {v c ) = ^a T cVU, &A Kv°~ 1 v~ a 5 pi 


(A4) 


(A5) 


Using (A2) to eliminate K, we obtain: 


5 = Bp 


(A 6 ) 


where the factor r/ is: 


V 


3 c(q) 1 L c [y c ) 

2 a T c UradZ /" -1 L s {y s ) 



l/(l+a) 


(A7) 
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• Equipartition condition. 

If we assume equipartition (Ub — U e ) between the magnetic energy density Ub and the 
energy density of relativistic electrons: 


/*7max 

U e = me 2 / ^( 7 ) 7 ^ 7 . 

J 7min 


we obtain the following relation between 5 and B: 


S = £ 5”( 3 +“)/p 


where £ is: 


(A8) 


(A9) 


£ = 


8nmc 2 L s (v s )f 1 1/p 


c(a)v s a V 

The factor / depends on the details of the electron distribution 


(A10) 


/ = 


n — 2 


( 7 2 T n - 7 2_n ) 
\ / mm / max / 


Equating Eq.(A9) and Eq. (A6) we can write the magnetic held as: 


- c\ p/(p + 3 +“) 

b= { i) 

JlJ 


and, inserting in Eq. (A6) 




p/(p+3+a) 


(All) 


(A12) 


(A13) 


The normalization K of the emitting electrons can be derived rearranging the expression 
of the synchrotron emission, Eq.(A2). 
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Figure Captions 


Fig. 1.— Profile of the X-ray (1 keV) and radio (4.9 GHz) flux along the jet of PKS 1136-135 
(data taken from Sambruna et al. 2005). 

Fig. 2. — Profiles of the relevant quantities (r, top panel, B and K, lower panel) for regions 
B-F as estimated from the radiative model. 

Fig. 3.— The parameter x (top), the comoving density (normalized at the initial value, 
middle) and the pressure (in units of 10 -12 dyn/cm -2 ) as a function of the Lorentz factor 
of the jet, calculated assuming the initial conditions inferred for knot C of the jet of 1136- 
135. Solid lines refer to the case of an initial pressure equal to that of magnetic held and 
relativistic electrons, while dashed lines report the results assuming a pressure ten times 
that the non-thermal one. The dotted line in the density panel visualizes the value of the 
density in the case of deceleration without entrainment. Points overimposed on the pressure 
curve indicate the value of the pressure provided by magnetic held and non-thermal electrons 
derived for the different portions of the how through the modelling of the emission (Fig.2). 
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